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ABSTRACT 

This work reports theoretical investigations of the surface motion 
in solids caused by reflection of an intense acoustic pulse at a pressure 
release surface done using the surface velocity reflection coefficient. 

A formula was derived for the surface velocity reflection coefficient 
which relates the surface particle velocity for an incident compression 
wave to the particle velocities of the incident and reflected waves. 

To obtain the surface profiles, two different assumptions were used. 
The first assumption was the conventional one. In this assumption, 
surface motion was neglected during the period that the pulse strikes 
the surface and it is called the fixed surface approximation. The 
second assumption is to consider the motion and distortions of the 
surface and is called the moving surface approximation. 

Four different materials were investigated and the results were 
compared for different materials and in the two different approximations. 
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I. INTRODUCTION 



Impulsive loads such as are produced by explosions and impacts may 
cause the failure of materials. An impulsive load is distinguishable 
from a conventional load primarily by its suddenness of application and 
its brevity of duration. Because of these loading characteristics, the 
reaction of a solid body to an impulsive load must be understood in 
terms of the propagation of waves in the body generated by the loading. 
If the loading is of sufficiently large amplitude, these waves produce 
extensive fracturing and large permanent distortions in the body upon 
which it acts. The phenomena associated with failure under impulsive 
loads are markedly different from those connected with failure under 
conventional loads. The reaction of materials to the extremely high, 
very transitory pressures imposed on them by impulsive loading is of 
great fundamental importance [Ref. 1]. 

The basic character of the fracture and plastic flow which are 
caused by impulsive loading is beginning to be understood [Ref. 1], 

A common fracture process results by the generation of shear waves 
whose amplitude exceeds the shear strength of the material in which 
the wave is propagating. Such shear waves are usually generated by 
the reflection of a strong compressive wave from boundaries. Another 
common fracture process results when the pressure release wave generated 
by the reflection of the compressive wave from a boundary has too great 
an amplitude. A third common failure process involves the superposition 
of shocks reflected from various boundaries. Thus, understanding the 
way in which intense waves damage solids requires an understanding of 
reflection of waves from boundaries. 
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In the limited theoretical discussion which has been given such 
processes to date [Ref. 5], the reflection process is discussed using 
expression obtained for small amplitude waves. Experiments show that 
this discussion is sufficient for qualitative, but not always quantita- 
tive, understanding of some, but not all, of the observed phenomena. 

The laws governing oblique reflection of elastic waves at a free 
surface are well known [Refs. 2, 3]. There is a particle motion 
associated with the reflection of these elastic waves that strike 
free surfaces obliquely. The object of this thesis is to study 
possible departures from the behavior predicted by the reflection 
formulae which result because of the finite amplitudes of waves 
involved by investigating the surface motion in solids caused by 
reflection of an intense acoustic pulse at a pressure release 
surface. 
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II. SURFACE MOTION DUE TO REFLECTION OF ELASTIC WAVES 



A. ELASTICITY 

Mathematical theories of elasticity are developed in the same 

general way. First, the notations of stress and strain are developed. 

Second, a stress-strain relationship between these quantities, or 

their derivatives, is assumed which idealizes the behavior of actual 

materials. Finally, using this relationship, equations of motion or 

equilibrium are set up which enable the state of stress or strain to 

be calculated when a body is subject to prescribed forces. 

In general, the stress on a surface element in a solid body 

does not act normally to the surface. It has components both normal to 

the plane and tangential to it. In Fig. 2-1, the coordinate system has 

three mutually perpendicular axes 0 , 0 , 0 . The stresses acting on 

x y z 

three planes normal to these axes which pass through a point P are con- 
sidered. There will be nine components of stress. These are denoted 
by a xx , Oyy, o zz , a xz , etc. The first letter in the suffix indicates 
the direction of the stress and the second one defines the plane 
in which stress is acting. If an infinitesimal rectangular paral lelepiped 
around P, with its faces normal to the axes, is considered, by taking 
moments, it may be seen that for equilibrium a = a , o = o , and 

Ajf J A At Cm A 

0 y z = a Z y. Thus there are only six independent components of stress. 
These completely define the stress at the point. 

The strain of the medium is the relative displacement of the material 
particles in the medium. A rigid translation or rotation of a medium 
is not considered a strain. In order to find the strain at the point, 
its change in position relative to adjacent points must be analyzed. 
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Figure 2.1. 



Stress components acting on an infinitesimal rectangular 
parallelepiped 
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Consider a point very close to P, which before displacement had 
coordinates (x + dx), (y + dy), (z + dz). Let the displacement which 
it has undergone have components (u + 6u, v + 6v, w + &d), where u, v 
and a) are the components of the displacement of point P. If 6x, 6y 
and 6z are sufficiently small. 



su - If {x * F 6y + « {z 



( 2 - 1 ) 



5v = “ + If + H sz 



S " ' 3x 6X + fy + fz 6Z 



The relative displacement of all surrounding points P may be found 

... ., , . .. . .... 3U 3U 3U 3 V 3 V 3V 3uj 3 w 3'w 

if the values of the nine quant,t,es w H 

at the point are known. These nine quantities are regrouped as follows: 
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The first three quantities, e , e , e correspond to the fractional 

yy zz 

expansion. The second three e , e and e correspond to the components 

y Z ZA Ajr 

of shear strain in the planes denoted by their suffixes. The last three 

w , u and to do not correspond to a deformation of the element around P, 
x y z 

but are the components of its rotation as a rigid body. 

Experimentally it is found that for most solids, the observed strains 
are proportional to the applied load, provided that the load does not 
exceed the elastic limit. This law is known as generalized form of 
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Hooke's law, and can be stated as: Each of the six components of stress 

is at any point a linear function of the six components of strain. This 
law may be written as 

°xx = c ll E xx + c 12 e yy + c 13 e zz + c 14 E yz + c 15 E zx + c 16 e xy 



(2-3) 



°yy = c 21 e xx + c 22 E yy + c 23 E zz + c 24 E yz + c 25 E zx + c 26 E xy 
°zz ’ c 31 e xx + c 32 E yy + c 33 E zz + c 34 E yz + c 35 E zx + c 36 E xy 
°yz = c 41 e xx + c 42 e yy + c 43 E zz + c 44 E yz + c 45 E zx + c 46 E xy 



= C cl e „ + C co e ww + C co e_, + C c/ 1 e w , + C cc e „ + C cc e. 



u zx 51xx 52 yy Wyz ^S^zx 56 xy 

°xy = c 61 G xx + c 62 G yy + c 63 e zz + c 64 G yz + zx + c 66 G xy 
The coefficients in these relations are known as the elastic constants 
of the material. 

In an isotropic solid, the values of the coefficients are independent 
of the set of rectangular axes chosen. For such a solid, there are only two 

independent elastic constants, x and p known as Lame's constants. All the 

elastic constants can be expressed in terms of the Lame's constants. 

Thus, the Lame's constants completely define the elastic behavior of an 

isotropic solid. For such a material, equation 2-3 may be written 



°xx = X4 + 2ue xx °yy = U * 2l,E yy °zz = U * 2ee zz 



(2-4) 



O = uc 

yz M yz 



a = uc 
zx M zx 



°xy pe xy 



Here & = e + e + e , and is called dilatation. It represents the 
aa yy zz 

change in volume of a unit cube. 

For convenience, four elastic constants are normally introduced. 
One is Young's modulus E, which is given by a / e so that: 

/V^ AA 
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(2-5) 



F _ p(3A + 2 u) 
A + y 



A second is Poisson's ratio v defined as -o/e so that: 

yy xx 



(2-6) 



V - 



2(A + y) 



A third is the bulk modulus k defined as the ratio between the applied 
pressure and the fractional change in volume when the solid is subjected 
to uniform hydrostatic compression so that: 



(2-7) 






Finally, the rigidity of shear modulus, which is identical with Lamd's 
constant y, is defined to be the ratio between the shear stress and the 
shear strain. 

B. EQUATION OF NOTION IN AN ELASTIC MEDIUM 

In order to obtain the equations of motion for an elastic medium, 
the variation in stress across a small parallelepiped with its sides 
parallel to a set of rectangular axes are considered (Fig. 2.2). The 
components of stress will vary across the faces; to obtain the force 
acting on each face, the value of stress at the center of each face 
times the area of the face has been taken. 

Six separate forces will act parallel to each axis. The resultant 
force acting in the x-di recti on is 

(°xx + “T; <5x)6y6Z - aSySz + 

XX 3x xx , 



+ ( CT y V + 3 °*y <Sy)6x6z - 6 5X6Z + 
^ 0 y X -7 
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Figure 2.2. Stresses acting on a small rectangular parallelepiped 
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By Newton's second law of motion, neglecting body forces such as 
gravity, this will be equal to 



(p Sx 6y6z) 



3 2 u 

3t 2 



where p is the density of the element and u is the displacement in the 
x-di recti on. Thus the equations of motion are 



(2-8) 
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C. ELASTIC WAVES 



Equations 2-8, 2-9 and 2-10 will hold whatever the stress-strain 
behavior of the medium. By substituting the equations 2-4 for the 
stress components in Equation 2-8, one obtains 



3 U _ 3 / , . 

p 77 ' « ( s 



2uc xx' 



+ ay <ue >w ) 



32 (l,E X2 ) 



By using the equations 2-2, one finds 

(2-11) p-Mf = (x + m)!^- + 

3t d 9X 
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Similarly 




Equations 2-11 through 2-13 are the equations" of motion of an 



isotropic, elastic solid in which body forces are absent. They predict 
the propagation of two types of waves through the medium. If both 
sides of the above equations are differentiated with respect to x, y and 
z respectively, they become 



This is the wave equation and shows that the dilatation A is propagated 
through the medium with velocity 

(2-15) c, = [(X + 2p)/p]^ 

On the other hand, if A is eliminated, it can be shown that 
■ Jl- 

(2-16) p — j- - yV a) 

3t 

where u> is the rotation about the x-axis. Similar equations may be 
obtained for u and w . Thus the rotation is propagated with velocity 



This shows that in the interior of any elastic solid waves may be 
propagated with two different velocities. These two types of waves 
should be termed dilatation waves and distortion waves. 



(2-14) 




= (X + 2y)v 2 A 



(2-17) c 2 = (v/p) h 
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Any plane wave propagated through on isotropic elastic medium 
must travel one or the other of the above velocities. The waves which 



propagate with velocity c-j are called longitudinal waves. The motion 
is along the direction of propagation. The other type waves, which 
propagate with velocity C£, are called transverse v/aves since the 
motion is parallel to the wave front. 

The equations 2-14 and 2-16 are of the form 



When the deformation is a function of only one coordinate, e.g. x, the 
equation becomes: 



where f and F are arbitrary functions depending on initial conditions, 
f corresponding to a plane wave travelling along the positive direction 
of the x-axis, and F to one in the direction opposite to this. 

D. REFLECTION AT PRESSURE RELEASE SURFACE 

Normal incidence . Any elastic wave will be reflected when it 
reaches a free surface of the material in which it is travelling. The 
simplest case occurs when the wave strikes the surface normally. In 
a longitudinal wave, since the stress normal to the surface at the 
surface must be zero, the reflected pulse must be opposite in sense 



(2-18) . ,2,2, 



at’ 



(2-19) 




The general solution of this is: 



i l> = f(x - ct) + F(x + ct) 
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to the incident pulse (compression reflected as tension and vice 
versa). 

Oblique incidence . When an elastic wave strikes a free surface 
obliquely, the situation is much more complicated. This complication 
arises from trie fact that the energy of incident wave is partitioned 
into two reflected waves instead of only one. When either a longitudinal 
or a transverse pulse strikes a surface obliquely, both a longitudinal 
and a transverse wave are reflected. (Figure 2.3). 

Consider first a longitudinal wave striking a surface obliquely. 

There will be a reflected longitudinal portion in which the angle of 
incidence a is equal to the angle of reflection. There will also be 
a transverse wave which is generated at the point of reflection. 

From Snell's law, the relationship between the angle of incidence a 
and the angle of reflection $ is 



5TnT ' cf = [20 - "I/O - **)]< 



(2-20) Sin a = . 

Cl n O r* 

'2 

From equation 2-20, it is seen that the angle of reflection 3 is the 
function of the incidence angle and the Poisson's ratio of the material. 
Figure 2.4 shows 3 as a function of a for several values of Poisson's 
ratio. 

This generation of a shear wave car, visualized in terms of 
particle velocities (Fig. 2.5). It can also be seen that a shear 
wave generated by the oblique reflection of a longitudinal wave is 
polarized in the plane of incident; that is, the particle velocity 
associated with this shear wave is normal to the direction of propaga- 
tion and directed in only one direction. The case of a compression 
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Figure 2.3 Oblique incidence of a plane longitudinal wave at a free 

surface 



20 



3 



in 

CD 

CD 

U 

CD 

CD 

-o 



cn 



+J 

v 

<D 



O 

CD 

r*“ 

CD 

5 




Angle of incidence a (degrees) 

Figure 2.4. Dependence of Angle of Reflection of Generated Shear Wave 
upon Poisson's Ratio and Angle of Incidence of Longitudinal 

Wave 



21 



Y 




Figure 2.5. Particle Notion Associated with an obliquely incident 

plane longitudinal wave 
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wave striking the free surface of a semi -infinite body will be 
considered. Let c be the instantaneous value of normal stress of the 
incident wave, a 1 the corresponding normal stress of the reflected 
longitudinal wave and t' the shear stress of the induced transverse 
wave. These stresses are all related by a coefficient of reflection R 
such that [Ref. 3] 

(2-21) a' = Ra 

(2-22) t' = C(R + 1) Cot2f3]a 

and the coefficient of reflection R is given by [Ref. 3]. 

2 

(2-2 3) R = tan g tanj?g - tan a 
tan 6 tan‘"26 + tan a 

R is a function of a and 6. Since R is a function of s» it is actually 
a function of the Poisson's ratio of the material. Figure 2.6 shows 
R as a function of a for several values of Poisson's ratios. 

A very similar series of statements can be made for the oblique 
reflection of a shear wave. By referring back to R as defined by 
equation 2-23. 

(2-24) t' = Rt 

(2-25) a' = [(R - l)tan 2 B]T 

where t, t 1 and a' are the instantaneous values of incident shear 
stress, generated shear stress and generated longitudinal stress, 
respectively, s now representing the angle of incidence. 

When an incident wave strikes a free surface obliquely, the 
surface will be displaced such that the angle of emergence a satisfies 
[Ref. 3, 4]. 
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Figure 2.6. 



Coefficient 

Incidence 



of Reflection as a Function of Angle of 
for Several Values of Poisson's Ratio 
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(2-26) 



o = 2g 



The angle of emergence is also related to rectilinear components 
of the particle velocity by 

(2-27) a = tan~^(v x /Vy) 

where v is the sum of the particle velocity components perpendicular 

A 

to the surface, v^ is the sum of the particle velocity vectors parallel 
to the surface as seen in Fig. 2.5. Again it can be seen that the 
angle of emergence is a function of both the angle of incidence and 
Poisson's ratio. 

E. SURFACE VELOCITY REFLECTION COEFFICIENT 

The surface particle velocity for an incident compression wave is 
related to the particle velocities of the incident and reflected waves. 
For a plane elastic longitudinal wave, the particle velocity at any 
time is given by the expression [Ref. 4] 

(2-28) = a/pc-j 

where a represents the normal stress associated with the wave. For 
a compressive wave the particle velocity is in the same direction as 
the velocity of wave propagation; for a tension wave the particle 
velocity is directed oppositely from that of propagation of the wave. 

A corresponding relationship exists for particle velocity within a 
transverse wave, namely, 

(2-29) = i/pC2 

where t is the shear stress associated with the transverse wave. 

From Fig. 2.5, the surface particle velocity is 
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(2-30) 



v p- v i* v t* v t 



By substituting equations 2-28 and 2-29 into the equation 2-30 

C -, 

(on, + o'rL + 

'1 1 1 v 2 



< 2 - 3 ') “p = Caft 1 + ° ,fi 2 + -ET T ‘ fi 3 ) 



where n-j is the unit vector in the direction parallel to the incident 
longitudinal wave, n 2 is the unit vector in the direction parallel 
to the reflected longitudinal wave and n^ is the unit vector in the 
direction perpendicular to the inflected transverse wave. By substitut- 
ing equations 2-21 and 2-22 into the equation 2-31 

Ci 



(2 ' 32) V^ [f 'l + fi 2 R+fl 3 4 (R+,)COt 



26 ] 



or 



(2-33) 



v * 



where 



(2-34) 



t = ^ Cn 1 + n 2 R + n 3 ^ (R + l)Cot 2 b] 



From Fig. 2.5 



‘1 



= 1 Cos a + 3 Sin a 



n n = 



n~ = 



-l Cos a + j Sin a 
i Sin 3+3 Cos 3 



By substituting these equations into equation 2-34, it can be shorn 



that 
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(2-35) 



t 
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pC-j 



j\*[(l - R) Cos a + 

+ 3CO + R) (Si n a + 



(1 + R) Cot 26 Sin 6] + 
c, 

-i- cot 26 Cos 6)] 
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2 




On the other hand 



(2-36) 



v = v (l Cos 26 + 3 Sin 26) 
P P 



By using equations 2-33, 2-35 and 2-36, it can be shown that 



(2-37) 



0 (1 + R) f Cot 26 Cos 6 + Sin 6 } and 

R = pc 2 1 STT26 ' 

V = Ra 

P 



where R is called the surface velocity reflection coefficient. From 
this formula, it can be seen that the surface velocity reflection 
coefficient is a function of the angle of incidence, the Poisson's 
ratio and the density of the material. 
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III. INVESTIGATION OF SURFACE MOTION 



A. NATURE OF INVESTIGATION 

The theories of oblique reflection of elastic waves at a free 
surface were outlined in the previous chapter. Such waves cause a 
particle motion on the surface. The motion of the surface will cause 
the reflection characteristics of an acoustic pulse to change as 
various portions of the pulse encounter to surface. The theoretical 
investigations of the surface motion in solid caused by reflection 
of an intense acoustic pulse at a pressure release surface were 
facilitated by using the surface velocity reflection coefficient R 
which was derived in equation 2-5. 

The intent was to illustrate the surface motion for different 
materials, for different loadings, for different angles of incidence. 

For this purpose the IBM-360 computer was used. The surface profiles 
were obtained by using two different approaches. Tne first approach 
was the conventional one which is called the fixed surface approxima- 
tion. This assumes that surface motion can be neglected during the 
period that the pulse strikes the surface. The second approach 
is consider the motion and distortions of the surface. This is 
called the moving surface approximation. The results of these two 
approximations were then compared. 

The computer programs used in this investigation are listed in 
Appendix A. The fixed surface approximation computes the surface 
distortion by simply integrating the surface velocity, obtained by 
using a constant surface velocity reflection coefficient, resulting 
from a pulse arriving at the x = 0 plane. The moving surface approxima- 
tion makes two major corrections. First, the surface velocity is 
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computed at each time from the displaced position of the surface, 
the displacement resulting from the previous portions of the pulse. 
Secondly, account was taken of changes in the angle of incidence 
caused by distortions of the surface from a plane. This was done 
by constructing the normal to an equation describing the surface 
obtained by passing a quadratic curve through the point in question 
and its neighboring points. Because of the time and position 
dependent changes in the pulse shape, numerical integration was 
used to convert the surface velocity into surface displacement. 

Throughout this investigation a saw-toothed shape stress pulse 
with 1 mm. pulse width has been used to illustrate the surface motion. 
The amplitude of the stress was changed from a value equal to the normal 
fracture stress to ten times this quantity. The normal fracture stress 
is the maximum tensile stress which the material will tolerate. A more 
exact definition would be the minimum dynamic tensional stress required 
to rupture the material. The angle of indidence was varied from 5 to 
85 degrees, using 5-degree increments. Four materials, aluminum, steel, 
plexiglass, and lead, were investigated. The characteristics of these 
materials are given in Table I. 

B. RESULTS 

Comparison of Approximations : As an illustration of the differences 

between the two approximations, the surface motion produced in aluminum 
and lead for a sawtooth compressive pulse of amplitude equal to the 
fracture stress incident at 30° computed in the two approximations is 
illustrated in Figs. 3.1 and 3.2. The maximum displacement produced 
by the pulse is the same in the two approximations. This result is 
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SOME PHYSICAL PARAMETERS FOR MATERIALS INVESTIGATED 
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iurface motions for aluminum under two different assumptions. 
\mplitude of stress is equal to normal fracture stress, a = 30°, 
time = 0.251 x 10-6 se c. 
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Figure 3.2. Surface motions for lead under two different assumptions. 

Amplitude of stress is equal to normal fracture stress, a = 30°, 

time = 0.741 x 10" 6 
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expected since the net displacement can be computed from the impulse 
contained within the incident pulse and the elastic moduli of the 
medium within which it propagates. The moving surface approximation 
predicts a more rapid displacement of the surface than the fixed 
surface approximation. This also is quite understandable. 

The Angle 4> : The change of the angle of incidence due to surface 
distortion is called During the motion of surface <t> at any given 
spot changes with time. To illustrate the change in <t> by the time at 
a given position, the curve is drawn for aluminum and given in Fig. 3.3. 
It shows that the angle <i> starts to increase after the pulse strikes 
to the surface, and reaches its maximum value, then decreases again. 

But the decrease is relatively slower than the increase. This is a 
result of the assumed sawtooth pulse. The angle 4> is also effected 
by the angle of incidence and by the amplitude of stress. The 
effect of the angle of incidence is illustrated in Fig. 3.4. 

Figure 3.4 shows the change in angle <j> by the angle of incidence for 
all materials investigated. The angle <j> changes linearly with stress 
amplitude because of using the Hooke's Law and the relatively small 
displacements. The change in the angle <f> with stress amplitude is 
given in Fig. 3.5. In this figure, aluminum is chosen as an example. 

In the moving surface assumption, the distortion of the surface 
from a plane effects the pulse shape and causes distortion on the pulse. 
Assuming plane wave propagation, there is no distortion at the front 
and tail end of the reflected pulse since it is reflected from a flat 
surface. But in between these two ends, the different portions of the 
pulse cause different distortions on the surface which result in 
different angle of correction <j>. Because of different <j>, the reflected 
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Figure 3.3. 



Change in angle <}>with time for aluminum, 
equals to normal fracture stress. 



Stress amplitude 
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Figure 3.4. The change in angle <j> with angle of incidence a for different 
materials. The amplitude of pulse is normal fracture stress. 
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(degree) 




Figure 3.5. The change in angle $ mx with stress amplitude for aluminum, 
a = 60°, t = 0.251 x 10-6 S ec. 
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angle is not the same as the instant angle of incidence for these 
portions of pulse. Thus the wave fronts for these portions of pulse 
are not in the same direction relative to the front and tail end of 
the pulse. Hence the reflected pulse under the moving surface 
assumption is not the same as the incident pulse. From the investiga- 
tion, it was seen that the pulse became far apart after the reflection 
from the surface. And this causes the distortion in the stress field 
which is not investigated in here. 

The moving surface assumption also effects the amplitude of the 
reflected pulse. The amplitude of the reflected stress pulse equals 
to the product of the coefficient of reflection and the amplitude of 
the incident stress pulse. But from Fig. 2.6, it can be seen that 
the coefficient of reflection is the function of the angle of incidence. 
Since the assumption causes a change in the angle of incidence, the 
amplitude of the reflected pulse is not the same as the amplitude of 
the reflected stress pulse under the fixed surface assumption condition. 
It is either weaker or stronger than it. It depends on the angle 
of incidence and the material used. To illustrate the effect of <f> 
on the amplitude of stress pulse, two numerical examples were done 
for aluminum. The coefficients of reflection were calculated by 
using the equation 2-23 for two different angles of incidence, 30 
and 40 degrees. The amplitude of stress was the normal fracture 
stress and the maximum angles of correction were used in both cases. 

For 30 degrees angle of indidence, by using fixed surface approximation, 
the coefficient of reflection equals -0.778. But in the same condition 
by using the moving surface approximation, it equals -0.790. For 
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45 degrees angle of incidence, the first approximation results in the 
coefficient of reflection equals to -0.579, the second approximation 
results in the coefficient of reflection equals to -0.590. 

Surface Displacement : Surface displacement is effected by the 

stress amplitude and by the angle of incidence. The surface displace- 
ment increases linearly with stress amplitude. However, it decreases 
by increasing the angle of incidence. The surface motions for two 
different materials, aluminum and lead, are illustrated in Figs. 3.6 
and 3.7. 
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Y(xlO m) 




Figure 3.6. Surface motions for aluminum. Stress amplitude equals to 
normal fracture stress and a = 30. 
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Figure 3.7. 



X(xl0" 8 m) 



Surface motions for lead 
normal fracture stress, 



. Stress amplitude equals to the 
a = 30°. 
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APPENDIX A 



THE COMPUTER PROGRAMS 
FOR 

SIMULATION OF SURFACE MOTIONS 
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